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PREFACE 


Under  the  sponsorship  of  the  Tactical  Technology  Office  of  the 
Defense  Advanced  Research  Projects  Agency,  The  Rand  Corporation  has 
been  engaged  in  analysis  for  and  development  of  hydrodynamic  design 
criteria  employing  various  concepts  of  boundary-layer  control,  includ- 
ing shaping,  suction,  and  heating. 

Stabilization  of  the  laminar  boundary  layer  over  a horizontal 
axisymmetric  body  by  heating  may  be  upset  by  the  buoyancy  force  induced 
by  temperature  differences.  Using  analytical  methods,  this  report  con- 
siders the  combined  effects  of  variable  viscosity  and  buoyancy  force 
over  a heated  cylinder,  as  a preliminary  step  toward  defining  a region 
where  heating  is  a favorable  method  of  stabilizing  the  boundary  layer. 

This  report  should  be  useful  to  hydrodynamicists , to  designers  of 
submersibles,  and  to  others  interested  in  applying  fluid  mechanics  to 
the  prediction  and  improvement  of  the  performance  of  underwater 
vehicles. 

Other,  related  Rand  publications  include: 


R 1752-ARPA/ONR,  Loui-Speed.  Boundary-Layer  Transition  Workshop 

William  S.  King,  June  1975. 

K-1789-ARPA,  Controlling  the  Separation  of  Laminar  Boundary 
Layers  in  Water:  Heating  and  Suction,  J.  Aroesty  and  S.’  A. 

Berger,  September  1975. 

R-] 863-ARPA,  The  Ljjects  oj  Wall  Temperature  and  Suction  on 
Lcvmnar  Boundary-Layer  Stability,  William  S.  King,  April  1976. 

R-1866-ARPA,  Hydrodynamic  Considerations  in  the  Design  of  Small 
Submersible  Vehicles  (U),  C.  Gazley,  Jr.,  J.  Aroesty,  W.  S. 
King,  and  E.  R.  Van  Driest,  April  1976  (Confidential). 

9 

R-1898-ARPA,  "e  " Stability  Theory  and.  Boundary -Layer  Transition, 

S.  A.  Berger  and  J.  Aroesty,  in  process. 

R- 190 7 -ARP A,  Buoyancy  Cross-Flow  Effects  on  the  Boundary  Layer 
oj  a Heated  Horizontal  Cylinder,  L.  S.  Yao  and  Ivan  Catton, 
April  1976. 

1. un-Shin  Yao  is  a member  of  Rand's  research  staff;  Ivan  Catton  is 
an  Associate  Professor  of  Engineering  at  the  University  of  California, 
Los  Angeles,  and  a consultant  to  I lie  Rand  Corporation. 
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SUMMARY 


The  stabilization  of  laminar  boundary  layers  under  realistic  cir- 
cumstances requires  attention  to  all  those  effects  that  could  tilt  the 
balance  towards  increased  growth  of  disturbances.  One  such  effect, 
occurring  in  heated  boundary  layers  in  water,  is  due  to  buoyancy- 
induced  cross  flow.  A small  cross  flow  is  induced  in  an  otherwise 
axially  symmetric  laminar  boundary  layer  when  a uniform  horizontal 
stream  flows  along  the  inner  or  outer  surface  of  a heated  horizontal 
cylinder.  The  magnitude  of  this  cross  flow  depends  on  the  ratio  of 
the  Grashof  number  to  the  square  of  the  Reynolds  number,  based  on  the 

radius  of  the  cylinder,  and  in  its  early  stages  grows  linearly  in  the 

2 

downstream  direction.  The  governing  parameter,  £ = Gr/Re  , is  inde- 
pendent of  viscosity  and  can  be  interpreted  intuitively  as  the  ratio 
of  the  heating-associated  buoyancy  force  to  the  inertial  force  of  the 
basic  flow. 

Earlier  analysis,  performed  for  a constant  property  flow,  indi- 
cated that  the  cross  flow  could  reduce  heat  transfer  and  shear  stress 
over  the  upper  half  of  the  cylinder  (or  the  pipe),  while  increasing  it 
over  the  bottom  half.  More  significantly,  the  cross-flow  profiles  were 
notoriously  unstable,  and  their  impact  on  hydrodynamic  stability  could 
be  important,  depending  on  the  magnitude  of  £. 

The  present  work  extends  the  analysis  in  our  earlier  study  to  in- 
clude the  viscosity  variation  appropriate  to  water  in  t lie  temperature 
range  40°F  (4°C)  to  100°F  (37.8°C).  The  earliest  effects  of  cross  flow 
are  studied  and  characterized  by  an  expansion  procedure  in  which  the 
basic  axially  symmetric  flow  determines  the  leading  term,  the  first 

corrections  to  axial  velocity  and  temperature  profiles  are  proportional 

2 

to  e(x/a)  cos  cf) , and  the  cross-flow  component  is  proportional  to 
e(x/a)  sin  (j). 

The  use  of  this  expansion  procedure  bypasses  the  need  for  numer- 
ical calculation  of  three-dimensional  boundary  layers,  provides  solu- 
tions which  are  valid  for  realistic  ranges  of  the  parameters,  and, 
because  of  their  explicit  dependence  on  the  cross-flow  parameter,  e, 
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can  be  later  studied  systematically  for  the  effect  on  cross-flow 
stability  of  increasing  the  magnitude  of  £.  However,  the  cross-flow 
effect  on  the  flow  transition  can  be  qualitatively  assessed  by  the 
variation  of  the  velocity  gradient  on  the  wall.  It  has  been  demon- 
strated by  Wazzan  and  Gazley  (see  Fig.  7)  that  the  critical  Reynolds 
number  can  be  correlated  approximately  with  the  velocity  gradient  on 
the  wall.  The  numerical  results  show  that  the  variable-viscosity 
effect  can  increase  the  velocity  gradient  and,  hence,  stabilize  the 
laminar  boundary  layer;  the  cross-flow  effect  will  decrease  the  veloc- 
ity gradient  and  destabilize  the  laminar  boundary  layer  over  the  upper 
half  of  the  cylinder  (or  the  pipe  flow).  Also,  the  boundary  beyond 
which  the  cross-flow  effect  can  overwhelm  the  variable-viscosity 
effect  has  been  determined. 

The  ratio  of  the  maximum  cross-flow  velocity  to  the  free-stream 
velocity,  u^,  is  about  3 x l(f3(AT/up  • x,  where  AT  is  in  °F,  u^  is 
in  ft/sec,  and  x is  in  ft. 

The  analysis,  which  may  easily  be  extended  to  larger  values  of 
£ (x/a) 2 , is  appropriate  for  estimating  cross-flow  effects  for  either 
external  or  internal  flows  and  will  be  reported  when  it  is  finished. 

It  is  thus  valid  for  describing  the  early  stages  of  buoyancy-induced 
cross  flow  within  a heated  tube,  where  experiments  are  being  conducted 
to  verify  the  theoretical  results. 

"A.  R.  Wazzan  and  C.  Gazley,  Jr.,  "A  Parametric  Study  of  Boundary- 
Layer  Stability  and  Transition  for  Falkner-Skan  Wedge  Flows,"  presenta- 
tion at  Lou-Speed  Bound<vv)-Uvjer  Transition  Workshop:  T1 , held  at  The 

Rand  Corporation,  Santa  Monica,  California,  September  13-13,  1976. 
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SYMBOLS 


a = radius  of  the  cylinder 
b,  c = constants,  Eq.  (5) 

f = nondimensional  stream  function,  Eqs.  (9),  (10),  and  (11) 

F = similarity  nondimensional  stream  function,  Eqs.  (14) 
and  (15) 


g = similarity  nondimensional  temperature,  Eqs.  (10)  and  (11) 

2 

g = gravitational  acceleration,  Eq . (2)  = 32.2  ft/sec 

G = similarity  nondimensional  temperature  independent  of 
Eqs.  (14)  and  (15) 

Gr  = Grashof  number,  Eq . (2) 

n = index,  Eq.  (17) 

N,  Nq,  = viscosity  ratios,  Eqs.  (2)  and  (6) 

Nu  = Nusselt  number,  Eq . (21) 

Pr  = Prandtl  number,  Eq . (2) 
r = nondimensional  coordinate  normal  to  the  wall 
r = radial  coordinate 
Re  = Reynolds  number,  Eq . (2) 

T = temperature 
u = axial  velocity 
v = circumferential  velocity 
w = radial  velocity 
x = axial  coordinate 

a = coefficient  for  variable  viscosity 

-4  o "1 

6 = thermal  expansion  coefficient  = 0.8  x 10  [ F] 

y = thermal  diffusivity 
e = Gr/Re^,  Eq.  (2) 

n = r//2x,  Blasius  similarity  variable,  Eq . (9) 

0 = nondimensional  temperature,  Eq . (2) 

V = kinematic  viscosity 

(jj  = circumferential  coordinate 

T = x-directional  shear  stress 
rx 

X = circumferential  shear  stress 
r<p 


Subscrip  ts 


0 = zeroth-order  solution 

00  = zeroth  order  of  zeroth-order  solution 

01  = first  order  of  zeroth-order  solution 

1 = first-order  solution 

10  = zeroth  order  of  first-order  solution 

11  = first  order  of  first-order  solution 

2 = first  order  cross-flow  solution 
w = surface 

00  = free  stream 
fc  = forced  convection 

Superscripts 

— = dimensional  quantities 
' = derivative  with  respect  to  r) 
x = local  quantities  at  location  x 
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I.  INTRODUCTION 


Laminar  flow  control,  as  the  means  to  achieve  low-drag  performance, 
requires  the  careful  manipulation  of  boundary-layer  velocity  profiles. 
Body  shaping,  pressure  gradient,  and  suction  are  classic  methods  for 
such  boundary-layer  control . ^ For  submersible  applications  the 
temperature-dependent  viscosity  of  water  provides  the  mechanism  by 
which  appropriate  velocity  profiles  can  be  maintained  by  wall  heating. 
These  velocity  profiles,  as  in  the  other  methods,  inhibit  the  amplifi- 
cation of  two-dimensional  Tollmein-Schlichting  instability  waves  and 
thus  retard  boundary-layer  transition  until  body  shape  and  adverse 
pressure  gradient  force  the  separation  of  the  laminar  boundary  layer. 
Suction  and  heating  have  similar  effects  on  boundary-layer  stability 

and  transition,  but  the  theoretical  potential  of  suction  for  even  f4Mfc 

(2) 

ther  drag  reduction  by  delaying  laminar  separation  is  greater. 

Historically,  each  method  of  boundary-layer  control  has  had  limi- 
tations, corresponding  to  the  growth  of  disturbances  associated  with 
free-stream  turbulence,  surface  roughness  or  waviness,  noise,  loss  of 
suction  effectiveness,  or  three-dimensionality.  Surface  heating  shares 
these  possible  limits  on  performance  and  has  some  limits  which  are 
unique  to  it  alone. 

Temperature  gradients  in  the  boundary  layer  of  a heated  horizontal 
body  result  in  buoyancy  forces  that  induce  small,  steady,  cross-flow 
velocities.  These  in  turn  could  distort  the  basic  laminar  axisymmetric 
velocity  profiles  into  three-dimensional  ones,  which  are  considerably 

(3) 

less  stable.  Yao  and  Catton  earlier  examined  the  development  of 
this  steady,  slightly  three-dimensional  laminar  boundary  layer  on  a 
heated  longitudinal  horizontal  cylinder  for  the  case  of  constant  fluid 
properties.  Their  results  indicate  that  the  cross  flow  can  enhance  the 
heat  transfer  and  stabilize  the  boundary  layer  over  the  lower  half  of 
the  cylinder,  but  it  degrades  heat  transfer  and  destabilizes  the  flow 
over  the  upper  half  of  the  cylinder.  Since  the  boundary  layer  is  very 
thin  compared  with  the  curvature  of  the  cylinder,  the  transverse  curva- 
ture effect  can  be  neglected.  Also,  the  results  can  be  applied  to  the 
entrance  flow  in  a circular  pipe. 


The  present  paper  extends  that  work  to  the  realistic  case  of 
water,  including  the  variable  viscosity  effect.  Similarly,  the  results 
can  be  applied  to  the  pipe  flow  near  the  entrance.  In  this  way  the 
domain  where  cross-fli  v effects  may  counterbalance  the  stabilizing 
variable  viscosity  effect  may  be  established  later  by  a more  detailed 
three-dimensional  stability  analysis. 

The  physical  model  chosen  is  a semi-infinite  cylinder  of  radius 
a (or  in  a horizontal  pipe  with  radius  a),  which  is  alignei  with  its 
axis  parallel  to  a uniform  flow  and  normal  to  the  direction  of  gravity. 
The  uniform  flow  is  assumed  to  have  a velocity  u^  and  temperature  T . 

The  surface  of  the  cylinder  is  heated  to  a constant  temperature 
T <T  > T ).  For  water  flow  in  the  range  of  40°F  (4.4°C)  and  100°F 
(37.8°C),  the  principal  departure  from  constant  property  flow  is  due 
to  viscosity  variation.  The  thermal  conductivity  and  specific  heat 
do  not  vary  appreciably  in  this  range,  and  the  most  important  phenomena, 
because  of  the  gradient  of  viscosity,  are  quite  well  represented  by 
the  model  in  which  k and  c^  are  constant,  and  the  variation  of  viscosity 
witli  temperature  is  preserved.  Since  the  density  variation  is  small, 
the  Boussinesq  approximation  is  adopted,  which  treats  the  density  as  a 
constant  in  the  equations  of  motion  except  for  the  buoyancy  force  term. 
The  buoyancy  cross  flow  is  small  in  the  leading-edge  region  of  the 

cylinder  and  can  be  treated  as  a second-order  effect.  Further  down- 

u 

stream,  a distance  of  order  a • Re/Gr2,  the  initially  small  buoyancy 
cross  flow  becomes  one  of  the  dominant  velocity  components  and  can  no 
longer  be  treated  as  a second-order  effect.  The  present  analysis  is 
valid  under  most  practical  circumstances,  particularly  when  realistic 
magnitudes  of  l//e  = Re/Gr  ~ = u^/VSg  A'fa  are  considered.  For  example, 
the  analysis  should  provide  an  accurate  approximation  to  the  cross 
flow  in  the  region  0 < x < a//e  = AT  • /a),  where  x is  in  ft, 

uoo  -*-s  in  ft/sec,  AT  is  in  deg  F,  and  a is  in  ft. 

fhe  ratio  of  the  cross-flow  component  to  the  free-stream  velocity 
is  Bg  ATx/u^  cos  if,  or  3 x 10  (AT/u^)  x x cos  ip,  in  the  same  units. 
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II.  ANALYSIS 


The  governing  equations  of  the  flow  with  fluid  of  variable  vis- 
cosity on  a heated,  longitudinal,  horizontal  cylinder  are  the  Boussinesq 
boundary-layer  equations.  In  cylindrical  coordinates,  as  shown  in 
Fig.  1,  they  are 
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Equations 


(1)  are  also  valid  for  the  entrance  flow  in  a heated  pipe. 


u 


CO 


Fig.  1 — Physical  model  and  coordinates 
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The  dimensionless  variables  that  are  introduced  to  nondimension- 
alize  Eqs.  (1)  are 


wkRe 


V = — , w = 

U U 


0 = 


T - T 

00 

T - T 
w ° 


x = - , r = 


u a 


Re  = 


C-r  = 


8ga3(T 


- T ) 


v 


Pr  = 


Y 


(the  velocities) 


(the  temperature) 


_ (r  a)^Re  (the  coordinates) 


(the  Reynolds  number)  (2) 


(the  Grashof  number) 


(the  Prandtl  number) 


£ 


N 


The  radius  of  the  cylinder,  a,  has  been  selected  as  the  longitudinal 
characteristic  length  in  order  to  study  the  detailed  development  of 
the  mechanism  of  the  induced  cross  flow  caused  by  heating  along  the 
upstream  part  of  the  horizontal  cylinder.  The  coordinate  normal  to 
the  wall  has  been  stretched  to  reflect  the  fact  that  the  thickness  of 
the  boundary  layer  in  the  region  of  x ~ 0(a)  is  inversely  proportional 
to  the  square  root  of  the  Reynolds  number. 

We  note  that  the  quantity  £ can  be  expressed  simply  as  Bg(Tw  * roo 
and  can  be  interpreted  as  the  ratio  of  the  potential  energy  increment 
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to  the  kinetic  energy  of  the  flow.  For  the  flow  close  to  the  leading 
edge  of  the  cylinder  (or  close  to  the  entrance  of  a pipe),  the  magni- 
tude of  e determines  the  importance  of  the  thermally  driven  cross  flow. 
As  we  move  further  downstream,  the  governing  parameter  will  be  the 
Grashof  number  instead  of  £.  There  the  cross  flow  can  be  neglected 
only  when  the  Grashof  number  is  small.  For  our  case,  where  the  Grashof 
number  is  not  small  but  Gr/Re  is  small,  the  cross-flow  effect  can  be 
treated  as  a secondary  effect  in  the  region  near  the  leading  edge  of 
the  cylinder  (or  near  the  entrance  of  the  pipe) . The  magnitude  of  the 
cross  flow  develops  and  eventually  becomes  one  of  the  dominant  phenomena 
which  can  trigger  the  flow  transition  and  the  flow  separation  further 
downstream.  it  is  worth  noting  that  the  value  of  £ also  determines  the 

size  of  the  early  region  where  the  cross  flow  is  a second-order  effect 

2 

when  Gr  > 1 and  Gr/Re  <1.  In  terms  of  the  dimensionless  parameters 
in  Eq.  (2);  F,qs.  (1)  become 


9u  9w  9v 

9x  9r  9c() 


0 


(3a) 


u 


9u 

9r 


- + f-  (»  |h> 

9x  9r  9r 


(3b) 


• 0 • COS  0 

3r  v^Re 


(3c) 


9v 

9x 


+ v 


9v 

9<j> 


+ w 


9v 

9r 


= E 


sin  0 • 0 


_ lE.  + 

90  9r 


« & 


(3d) 


30.  lix  90  1 920 

9x  V 90  W 9r  Pr  ^2 


(3e) 


after  neglecting  smaller-order  terms.  The  terms  that  represent  the 
transverse  curvature  effect  have  been  neglected  simply  because  the 
boundary  layer  is  thin  compared  with  the  radius  of  the  cylinder  when 
the  Reynolds  number  is  not  small.  Equation  (3c)  indicates  that  the 
pressure  gradient  normal  to  the  wall  is  negligible  to  the  lowest  order, 
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and  the  pressure  gradients  parallel  to  the  wall  can  be  evaluated  from 
the  inviscid  solution  at  the  edge  of  the  boundary  layer.  For  a uni- 
form free  stream,  these  pressure  gradients  vanish. 

It  is  necessary  to  specify  the  viscosity  ratio,  N,  before  the  * 
solutions  of  Eqs . (3)  can  be  presented.  It  has  been  shown  by  Gazley 
that  a good  approximation  for  the  kinematic  viscosity  of  water  in  the 
range  of  temperature  between  40°F  (4.4°C)  and  100  F (37.8  C)  is 

10  5 (4) 

V " -0.0807  + 0.0126T 

where  T is  in  °F  and  V is  in  ft2/sec.  This  suggests  a viscosity- 
temperature  model  of  the  form 


1 

y b +'cT 


(5) 


For  such  a model, 


i = ——  = i + cx  at  • 0 

N y 


where  AT  = T - T^,  and  a = c/ (b  + cTj . For  water  in  the  temperature 
range  between  40°F  (4.4°C)  and  100°F  (37.8°C),  a « 0.0151  (°F)  . It 

has  also  been  shown  by  Gazley  that  the  Prandtl  number  of  water  in  the 
range  of  temperature  40°F  (4.4°C)  and  100°F  (37.8°C)  can  be  approxi- 
mated by 


Pr 


455 


(7) 


where  T is  in  °F.  The  illustrative  calculations  described  below  were 

CO 

performed  using  a value  of  Pr  = 8. 


"Personal  communication  from  Carl  Gazley,  Jr.,  The  Rand  Corpora 
tion,  1976. 
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The  solution  of  Eqs.  (3)  with  the  viscosity-temperature  model, 
Eq . (6),  can  be  expanded  into  a series  in  £,  if  £ is  small,  so  that 


u = Uq  + £u^  + 


v = ev  ^ + ... 


w = Wq  + £w^  + 


0 = 0Q  + £01  + 


N = Nq  + £N1  + 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 


where 


N 


0 1 + a AT  0 


0 


(8f ) 


-a  AT 


(1  + a AT  0O) 


2 * 6!  = A(V  ‘ 91 


(8g) 


Substitution  of  the  expansion  given  by  Eqs.  (8)  into  Eqs.  (3), 
and  the  collection  of  terms  of  equal  order  will  result  in  the  pertur- 
bation equations.  The  perturbation  equations  of  lowest  order  are 


(eU):  W+VJ-O 


9”  + Er  fo0’  = 0 


(9a) 

(9b) 


where  the  prime  denotes  a derivative  with 
being  the  Blasius  similarity  variable.  I 
fined  by 


1 1 1 i to  n,  with  H = r //2x 
•u  is  ii  function  f^  is  de- 


U0  = fo 


(9c) 


WQ  " ^ (nf0  V 


(9d) 
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The  second-order  perturbation  equations 


are 


<e  ):  <»„£")'  + f0f”  - 4£’f;  + 5fjf1  + f" 


3f, 


= (A  • g • f|j)'  (10a) 


(N0f 2*)  ’ + f0f2  ~ 2f0f2  + sin  <J)  • 0Q  = 


(10b) 


k s" + V - + 5%h  + % 


= 0 


(10c) 


The  stream  functions  and  f2  and  the  temperature  g are  defined  by 


Uj_  = (2k)  2f  ’ (n,cj)) 


V1  = (2x^f2(n*(()) 


w,  = 


(2x)3/2  ( nf| 


3f, 


5f i 

i d<p 


o = (2x)“g(n,(j)) 


The  boundary  conditions  associated  with  Eqs.  (9) 


are 


f0  - fo  = 0 and  8o  = 1 at  n = o 


fo+  1 


and  0Q  -*•  0 as  n 


and  the  ones  associated  with  Eqs.  (10)  are 


9f 


fl  [2  ~ fl  ~ d^T  “ 0 » 8 = 0 at  n = 0 


fl  f2  0 ’ 


g = 0 as  r| 


f;  = 0 


at  <j>  = 0,tt 


(11) 


(12) 


(13) 


Equations  (10)  are  separable  in 
ables  of  the  form 


terms  of  t)  and  (})  with  dependent  vari- 


f (n,<j>)  = F1ln)  cos  <j> 

f2(n,4>)  = F2(n)  sin  ♦ 
g(n,<t>)  = c(n)  cos  <J> 


(14) 


Substituting  Eqs.  (14)  into  Eqs.  (10)  gives 


<Vl>'  + Vi  - 4foFi  - 5fSFi  + fSF2  ' (A  ' G ' f0> 


(N0F"C  + £0F”  - 2£’F'  = -0„ 


L c..  + foG.  - - -ej(5Fl  + f2) 


(15a) 

(15b) 

(15c) 


The  associated  boundary  conditions  for  Eqs.  (10)  are 


\ 

F (0)  = F|(0)  = F|  (°°)  = 0 
F (0)  = F^(0)  = F 2 (°°)  = 0 
G(0)  = C.  (°°)  = 0 


(16) 


The  solutions  depend  on  temperature  only  through  the 
are  integrated  numerically  for  a • AT  = 0,  .5,  1.0. 
this  corresponds  to  AT  = 0,  33  , 66  . 


factor  a AT  and 
If  T = 57°F, 

OO 


The  form  of  the  solution  suggests 
expansion  procedure  along  the  lines 


that  a continuation  of 


this 


-lo- 


ll = un (n » 4>) (ex2)n 


= ^ 0 (o,4>)(ex2)n 

t—i  n 


(17) 


= ex  vo(n,0) (ex2)n 


would  be  the  appropriate  continuation  scheme  further  downstream.  Con- 

2 

vergence  of  this  expansion  presumably  requires  that  ex  be  less  than 

unity,  and  the  accuracy  of  the  two-term  expansion  for  u and  0 and  the 

2 * 

one-term  expansion  for  v requires  that  ex  be  small. 


Preliminary  numerical  computations  at  Rand  of  the  three- 
dimensional  momentum  integral  boundary-layer  equation  suggest  that 
the  solution  given  by  Eq.  (17)  is  valid  for  0 < Cx^  < .2.  However, 
the  limits  of  validity  of  the  analysis  and  regions  where  the  cross 
flow  may  balance  or  even  dominate  the  variable-viscosity  effect  can 
be  better  determined  after  more  computations  of  the  downstream  three- 
dimensional  boundary-layer  flow  are  completed. 
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TTT.  NUMERICAL  RESULTS 


The  physical 
pressed  in  terms 
They  are 


quantities  of  velocity  and  temperature  can  be  ex- 
of  the  functions  defined  in  Eqs . (8),  (9),  and  (11) 


u 


f^  + £(2x)2F^ 


COS  <})+..• 


(18a) 


v = e(2x)  • * sin  4>  + • • • 


(18b) 


= -L-  (nf>  - f ) + £(2x)3/2(nF’  - 5F1  - F2) 

j U vJ 

vlx 


COS  <j>  + • • • 


(18c) 


0 


0Q  + e(2x)2  • C 


cos  4>  + • • • 


(I8d) 


The  local  shear  stress  at  the  cylinder 
axial  and  cross-flow  velocity  components: 


can  be  computed  from  the 


T 

rx 


and 


Introducing  the  series  e 


xpansion  (18),  this  can  be  written  as 


T 7 F'l,(0) 

ri_,1  + E(2x)  -cos**.. 

'‘rx'fc  0 


(19) 


for  the  axial  shear  and 


X ~ e/2x  F"  (0)  sin 
ref)  z 


(20) 


for  the  circumf erentia 
forced  convection  with  variable  viscosity. 


1 component , where  Ec  denotes  the  voice  for 
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The  local  Nusselt  number  can  be  written  as 


Nu 


Nu. 

fc 


1 + e(2x)‘ 


G’  (0) 

0J(°) 


COS  (J>  + . . . 


(21) 


Table  1 shows  the  values  of  the  quantities  F^'(0)/fp(0),  G ' (0) /6q(0) , 
and  F(j(0).  It  is  worth  noting  that  the  relative  modifications  of  the 
axial  stress  and  heat  transfer  caused  by  cross  flow  are  both  small 
(<9  percent)  and  insensitive  to  wall  temperature  in  the  range 
0 <;  a AT  £ 1,  or  up  to  AT  ~ 66°F.  The  circumferential  velocity  gra- 
dient, on  the  other  hand,  shows  a considerable  variation  with  wall 
temperature,  increasing  50  percent  in  the  range  between  a = 0 and 
a AT  = 1.  Interestingly,  the  increase  in  the  coefficient  F'^CO)  is 
balanced  by  the  decrease  in  the  wall  viscosity  caused  by  heating,  and 
the  quantity  T , the  product  of  the  velocity  gradient  and  the  vis- 
cosity at  the  wall,  actually  decreases  by  25  percent  in  the  range  be- 
tween a = 0 and  a AT  = 1. 


Table  1 

ENHANCEMENT  PARAMETERS  FOR  HEAT  TRANSFER 
AND  WALL  SHEAR  STRESS 


c' (o)/e^(o) 

F'i  (0)  / f q(0) 

F”(0) 

a = 0 

0.06981 

0.04905 

0.41499 

a AT  = 0.5 

0.06559 

0.04897 

0.53625 

a AT  = 1.0 

0.06296 

-i 

0.04899 

1 

0.63888 

From  Figs.  2,  3,  and  4 we  observe  that  the  reduction  of  water 
viscosity  caused  by  heating  modifies  the  axial  flow:  The  axial  flow 

is  accelerated  by  the  induced  cross  flow  over  the  lower  half  of  the 
cylinder  (-tt/2  < (J>  < tt/2)  and  decelerated  over  the  upper  half 
(tt/2  < (j>  < 3 tt/2)  . The  cross  flow  is  accelerated  by  buoyancy  from  the 
lower  stagnation  point  (tf)  = 0)  to  its  maximum  value  at  <J>  = tt/2  and 
then  decelerated  to  its  upper  stagnation  point  ((J)  = tt)  . The  decrease 


\ 
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of  viscosity  by  heating  enhances  the  magnitude  of  the  cross  flow  (1 
and  shifts  the  point  of  maximum  velocity  to  the  wall  as  shown  in  Fig.  4. 
Therefore,  the  cross-flow  velocity  gradient  at  the  wall  magnified 
appreciably  by  a slight  variation  of  the  fluid  viscosity.  The  induced 
axial  velocity  becomes  more  full  in  the  wall  region,  as  wall  temper- 
ature increases.  The  key  differences  between  the  cross-flow  velocity 
F£  and  the  co-induced  axial  velocity  Fj  is  that  the  maximum  value  of 
F'  increases  from  0.118  to  0.107  as  a AT  increases  from  0 to  1,  while 
the  maximum  value  of  Fj  decreases  from  .021  to  .015  over  the  same  range. 
Since  the  maximum  value  of  F|  is  so  small,  this  suggests  that  the  ef 
feet  of  buoyancy  on  the  boundary-layer  displacement  and  momentum  thick- 
nesses, 5*  and  0,  would  be  negligible  in  this  region. 


1 


IV . BUOYANCY  VERSUS  VARIABLE  VISCOSITY 


The  primary  importance  of  cross  flow  is  the  stability  of  the 
resulting  three-dimensional  boundary  layer,  and  this  will  be  consid- 
ered in  subsequent  research.  However,  it  is  still  of  some  interest 
to  compare  the  magnitudes  of  the  buoyancy  effect  with  the  variable 
viscosity  effect,  in  order  to  further  delineate  the  region  where 
buoyancy  effects  on  shear  stress,  heat  transfer,  and  separation  may 
safely  be  neglected. 

While  the  solutions  of  Eqs.  (9)  and  (15)  correspond  to  the  most 
general  case  of  arbitrarily  large  viscosity  variation,  explicit  esti- 
mates of  the  region  where  buoyancy  effects  on  heat  transfer  and  axial 
velocity  profiles  start  to  be  important  can  be  obtained  by  a further 
approximation. 

If  the  viscosity  variation  is  small,  because  of  small  values  of 
a AT,  the  effect  of  variable  viscosity  can  be  separated  out  from  the 
first-order  Eqs.  (9a)  and  (9b)  by  a linearization  procedure.  When 
a = 0,  the  solution  of  Eq.  (9a)  corresponds  to  the  Blasius  solution. 
The  expansion  procedure  becomes 


fo  = foo  + (a  AT)  * f 01  + 


0o  = eoo+(a  AT)  *eoi  + -’ 


(22) 


Substitution  of  Eqs.  (18)  into  Eqs.  (9a)  and  (9b),  or  collecting  terms 
of  equal  order,  gives  the  equations  of  the  zeroth  order 

(23a) 
(23b) 

and  the  equations  of  the  first  order 


f!"  + f" 

r00  00  00 


= 0 


— 0" 
Pr  00 


+ fooeoo 


= 0 


{*} 


• . - ,i  j'-  s 1 x-r  - -■  j.2,  » . ■ - . i »;  va. 
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f "'  + Cf  f"  + f"  f ) = (0  f"  ) ' 

01  U00  01  oo  or  ^ oo  oo' 


Pr  G01  + f00801  + f01600  ° 


The  associated  boundary  conditions 


foo(0) 

= foo(0) 

ii 

T 

O 

- O 

<4-1 

o 

II 

O 

O 

O 

CD 

= 1 , 

eooM  ■ 

foi<0) 

= f01(0) 

- V">  ■ 0 

O 

T— 1 

o 

CD 

= 0O1(«) 

= 0 

(24a) 

(24b) 


(25) 


The  solution  of  Eq.  (23a)  is  the  Blasius  solution;  the  solution  of 
Eq.  (23b)  is  the  forced-convection  energy  equation.  The  functions  fQ^ 
and  0Q1  represent  the  effect  of  the  variable  viscosity  on  the  forced 
convection.  The  numerical  values  of  and  fQ1  can  be  obtained  by 

integrating  Eq.  (24)  numerically.  Strictly  speaking,  the  above  ap- 
proach is  valid  only  for  small  values  of  a AT.  However,  comparison 
with  the  solution  of  Eq . (9)  shows  that  this  linearization  is  accurate 
in  the  range  0 < a AT  < 0.5. 

The  decrease  of  the  water  viscosity  caused  by  heating  will  in- 
crease the  axial  velocity  gradient  on  the  wall  and  consequently  will 

(23) 

tend  to  stabilize  the  boundary  layer.  ’ For  a horizontal  cylinder, 
the  cross  flow  will  be  induced  by  wall  heating  as  a by— product,  which 
will  decrease  the  axial  velocity  gradient  on  the  wall.  In  general, 
the  stability  of  a boundary  layer  is  related  to  this  velocity  gradient 

(4) 

on  the  wall  as  has  been  demonstrated  by  Wazzan  and  Gazley,  Jr.  and 
shown  in  Fig.  6.  Based  on  this  observation,  the  cross  flow  seems  to 
introduce  a destabilizing  el  feet  on  the  boundary  layer.  Also,  Llie 
inflection  point  of  the  cross-flow  velocity  profile  could  be  unstable. 
The  conditions  leading  to  flow  transition  can  only  be  deiined  by  ex- 
tensive stability  analysis.  Before  we  can  exercise  such  a complicated 
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Fig.  6 — Critical  Reynolds  number  versus  wall  velocity  gradient 
(/Sis  related  to  the  wedge  angle) 


three-dimensional  stability  analysis,  we  will  apply  a simple  criterion, 
the  variation  of  axia]  velocity  gradient  on  the  wall,  to  study  the 
tendency  of  stabilizing  or  destabilizing  the  boundary  layer  by  heating. 

The  purpose  of  the  study  i • rhis  section  is  to  provide  some  in- 
sight into  the  balancing  effects  of  variable  viscosity  and  cross  flow. 
Thus  the  discussion  is  limited  to  the  case  of  small  values  of  a AT  in 
order  to  supply  a simple  criterion  to  define  a region  and  its  functional 
dependence  that  heating  is  a practical  way  to  stabilize  the  boundary 
layer . 

The  double  expansion  of  velocity  and  temperature  in  terms  of  e 
and  a AT  can  be  written: 


I 


IrtW 


19- 


u = fJQ  + (a  AT) f + e(2x)2[F|Q  .+  (a  AT)F^J  cos  <f>  + • • • 
0 = e00  + (a  AT) 0Q1  + e (2x) 2 [Gq  + (a  AT)G^ ] cos  <t>  + • • • 


(26a) 

(26b) 


where 


Fi  = Fio  + (aAT)Fn  + ••• 

G = GQ  + (a  AT)G1  + . . . 


(27) 


and  F ^ and  G^  are  the  first-order-solution,  cross  flow-induced  quanti 

ties,  with  a constant  property  assumption. 

The  local  shear  stress,  Eq.  (19),  can  be  further  expanded.  It 

becomes 


r x 


(Wfc 


f"(°)  2 F10(0) 

= 1 + (a  AT)  TTTTTfT  + e(2x)  — * cos  <J>  + 


foo(0) 


fSo<0) 


(28) 


where  fc  denotes  the  case  of  forced  convection  with  constant  viscosity. 
Since  the  wall  viscosities  are  the  same,  Eq . (28)  can  be  viewed  as  the 
ratio  of  the  velocity  gradients  at  the  wall.  Values  of  ^(0) ^f00^0) 
and  F'^q (0) / f^n (0)  are  given  in  Table  2,  and  ^(0) /fo0(0)  ls  about 


L00 


twenty  times  bigger  than  FjQ(0) /f^(0) . This  indicates  that  the  variable- 
viscosity  effect  overwhelms  the  cross-flow  effect  over  the  most  forward 
part  of  the  cylinder.  However,  the  magnitude  of  the  cross-flow  effect 


Table  2 

RATIOS  OF  HEAT  FLUX  AND  SHEAR  STRESS 


f^(0)/E”0(°) 

>V0)/E5o(0) 

9oi(0)/eoo(0) 

c'0m/%0m 

0.89169 

0.04905 

0.17791 

0.06981 

ar.  i*. >■  *£ ’t.(v . : i- Afi. 
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increases  as  ~ x2  downstream.  The  variable-viscosity  effect  is  even- 
tually balanced  by  the  cross-flow  effect  far  away  from  the  leading 
edge.  Equation  (28)  clearly  points  out  the  location  where  the  effect 
of  stabilizing  the  boundary  layer  by  heating  will  be  balanced  by  the 
destabilizing  cross-flow  effect  over  the  upper  half  of  the  cylinder 
(tt/2  < <j>  < 3tt/2).  This  location  can  be  found  from  Eq.  (28)  by  equat- 
ing the  second  and  the  third  terms  on  the  right-hand  side,  which  yields 


x = 


(a  AT)  f01(0) 

h 

L £ Fio(0)_ 

\ 

n 

(29) 


COS 


Equation  (29)  is  plotted  in  Fig.  7,  which  shows  that  the  region  of  sta- 
bilizing the  boundary  layer  by  heating  shrinks  as  x increases.  it 
should  be  noted  that  the  expansions  (23)  are  valid  only  up  to  x ~ L(l//e) 
and  then  the  originally  small  cross  flow  becomes  one  of  the  dominant 
velocity  components  beyond  x ~ 0(l//e)  when  Cr  is  not  small.  Equation 
(29),  however,  indicates  that  the  merit  of  the  heating  to  stabilize  the 
boundary  layer  starts  to  fade  away  at  x ~ 0(a  AT/e)  “ and  disappears 
before  the  magnitude  of  the  cross  flow  becomes  appreciable,  if  0.  Al  is 

small. 

The  explicit  expression  of  Eq . (29)  in  terms  of  physical  quantities 
can  be  written  as 


x ~ _ ~ 
\ a 


a 


aBg  cos  <j> 


(30) 


Equation  (30)  reveals  that  the  size  of  the  stabilized  region  by  heating 
is  proportional  to  the  free  stream  speed,  is  inversely  proportional  to 
the  square  root  of  the  cylinder  diameter,  and  is  independent  on  AT. 

For  a cylinder  of  1 ft  in  radius  submerged  in  60°F  water,  x is  about 
1.5  at  uro  = 3 ft/sec  and  150  at  u^  = 30  ft/sec  for  4>  = tt.  This  shows 
that  stabilizing  the  boundary  layer  by  heating  a small  body  is  highly 
effective.  However,  extrapolating  the  test  data  of  a small et  body  to 
a larger  size  appears  to  be  dubious.  Furthermore,  the  existence  ot  the 
inflection  point  in  the  cross-flow  velocity  could  cause  earlier  separa- 
tion and  shift  the  transition  point  upstream. 
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\ Boundary  of  zero  first-order  stress 
\ 


2 3 4 5 6 7 


X ‘ V a AT 

Boundary  where  the  variable- viscosity  effect  balances 
the  cross-flow  effect 
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